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Abstract- Analytical and numerical models are developed for the microwave 
small-signal performance, such as transconductance, gate-to-source capaci- 
tance, current gain cut-off frequency and the optimum cut-off frequency of the 
AlGaAs/GaAs High Electron Mobility Transistor (HEMT), in both normal 
and compressed transconductance regions. The validated I-V characteristics 
and the small-signal performances of four HEMTs are presented. 

Nomenclature 

L : Gate length. 

Z : Gate width. 

pi : Low field mobility of AlGaAs layer. 

P 2 '• Low field mobility of two-dimensional electron gas. 

d : Thickness of AlGaAs layer. 

di : Thickness of undoped AlGaAs layer. 

w : Width of undepleted region in AlGaAs layer. 

Nd : Doping concentration of AlGaAs layer. 

n M : Sheet concentration of two-dimensional electron gas. 

n to : Equilibrium Sheet concentration of two-dimensional electron gas. 

ti : Permittivity of AlGaAs. 

E c i : Saturation electric field of AlGaAs. 

E c 2 : Saturation electric field of two-dimensional electron gas. 
v, : Saturation velocity of two-dimensional electron gas. 

/3 : Charge control coefficient. 

6 : Effective width of conduction channel. 

V t ho • Threshold voltage for two-dimensional electron gas. 

Vbi : Built-in voltage of Schottky gate on AlGaAs layer. 

V p : Effective pinch-off voltage of AlGaAs layer. 

1 Introduction 

High frequency solid state technology has been moving towards the use of the high electron 
mobility transistors in microwave and in high speed digital circuits because of its high 
frequency operation and of its tolerance to many forms of radiation. Several workers 
have been studying the GaAs HEMTs both theoretically and experimentally since its first 
introduction in 1980 [1]. Over the past years, analytical, numerical and/or computer-aided 
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Figure 1: Schematic Diagram of a Uniformly Doped AlGaAs/GaAs HEMT. 


models have been reported by many authors. But, because of the complexity in structure 
of this device, VLSI circuit designers demand a more accurate and compact model for theif 
design. 

Among other workers, C.Z. Cil and S. Tansal [2], in i985, proposed an analytical model 
which used the simple Trofimenkoff-type velocity-field linear relation [3]. Their modeled 
results agree very well with the experimental data. However, their model is good only fof 
the linear normal transconductance region; it does not cover the current saturation region 
and also the parasitic conduction in the AlGaAs layer. But the computer-aided design 
and simulation of the HEMT circuits demand a complete and more accurate model. In 
1986, G.W. Wang and W.H. Ku [4] developed a compact but complete analytical model 
which covers the whole operation range of the dc characteristics. This model calculates 
the I-V characteristics of four different HEMTs and compares the modeled results with the 
experimental data. We have chosen their model as the basis for this work and from this 
model we have developed analytical and numerical models to calculate the small-signal 
performances, such as transconductance, g m , gate-to-source capacitance, C gt) current gain 
cut-off frequency, fa, and the optimum value of the cut-off frequency, fr(opt) before current 
saturation occurs. 


2 DC Model 

The basic structure of a HEMT device is significantly different from a conventional field 
effect transistor. A cross sectional view of a uniformly doped AlGaAs/GaAs HEMT device 
is shown in Figure 1 . 

At low gate voltage, it has only one current conduction channel but at high gate voltage, 
it has two conduction channels: one is the two-dimensional electron gas (2-DEG) in the 
interface between AlGaAs and GaAs and the other is the parasitic conduction through the 
undepleted n"*"-AlGaAs layer. If the AlGaAs layer is not fully depleted by the Schottky 
gate and the heterojunction, then the free carriers under the gate are the two-dimensional 
electrons and the free electrons in the AlGaAs layer. The width of the undepleted AlGaAs 
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region can be approximated by [5] 

(1) 

By setting w = 0, the AlGaAs layer is completely depleted, one can obtain the critical 
value of the gate voltage, Vg as 

V ' = Vc{w = 0) = Vbi - q ^(d-d i -’±y ( 2 ) 

The V G < V c defines the normal transconductance region where only the 2-DEG is 
the current conduction channel and the Vg > V c defines the compressed transconductance 
region where both the 2-DEG and the undepleted AlGaAs layer are the current conduction 
channels. 

According to the charge control model [6], the sheet charge density of the 2-DEG can 
be approximated as a linear function of gate voltage and channel voltage which is given by 

n.(x) = (3{V a - V(x) - V th ) (3) 


where x is in the direction along the heterojunction. 

In this dc model, for mathematical simplicity, the Trofimenkoff-type [3] electron velocity- 
field relation has been used for both the 2-DEG channel and the AlGaAs parasitic con- 
duction channel. The linear electron velocity-field can be related as 


•u(x) 


pEjx) 


( 4 ) 


Here, E(x) is the electric field in the 2-DEG channel or in the undepleted AlGaAs layer 
and E c is the field at which the velocity of electrons reach the maximum value (saturation 
velocity). 

Using the charge control concept and the velocity-field relationship described above, 
the current conducting through the 2-DEG channel can be determined by 


h-DEG = Zqn,(x)v(x) 


( 5 ) 


Similarly, the current through the undepleted AlGaAs layer can be determined by 

iAlGaA s = ZqN d w(x)v(x) (6) 

Here, for simplicity, full ionization of the donor atoms has been assumed for the current 
through the AlGaAs layer. 


(A) I-V Equations in the Normal Transconductance Region 

When the gate voltage is low, i.e. Vq < V c , the normal transconductance region is formed. 
This region is then divided into the linear (Vo < V sa t) and the saturation ( Vb > V sa t) 
regions. The current-voltage relationship in two different regions can be derived as follows: 
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Figure 2: Schematic Diagram Showing Current Saturation in the 2-DEG Channel. 


(i) Normal Linear Region (V# < U, ot ) 

Introducing equations (3) and (4) in equation (§), and integrating from source to drain 
along the 2-DEG channel, the current through the channel is 

My a - v* - %)v D 

lD ~ — r +> < 7 > 

where A = an d B — LE c2 are the model parameters; V G and Vb are the internal 

gate and drain voltages. 

(ii) Normal Saturation Region (Vjp > V aat ) 

The velocity-field relation (equation (4)) allows the velocity to saturate when the electric 
field approaches infinity. But physically it is impossible; so, the model assumed that the 
velocity saturation occurs when E > E c . Therefore, from equation (4), the saturation 
velocity at E = E e is v, = 

When the drain voltage, Vjp becomes greater than the saturation voltage, V, at the 
situation becomes like Figure 2. At x = Z c , electric field exceeds saturation field, E e , 
and the electron velocity saturates; and after this the electrons move with this constant 
saturation velocity. Then, using V = V sat and ^ = E c2 at x = L c , equation (7) can be 
written as 

r _ ZqMv G - Vtn - 

i. + fc () 

Also from equations (3) and (5), the current in the saturation region can be written as 


I D = Zq[3(V G - V th - V, al )v, = 


- V th - V, at )E c2 


Now, using the current continuity condition, equations (8) and (9) can be combined to 
obtain 

(l-g,)B(V 0 -V, t ) 

(1 - K,)B + (V a - V lh ) ll0) 


ilium ii in 
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where K\ = - . Generally, L c and V, at can be determined by solving a two-dimensional 

P oisson’s equation which has the form in the velocity saturation region 


d 2 V 

dx 2 


« aI D 


( 11 ) 


where a = • Here, S is the effective width of the conduction channel which is assumed 

to be invariant to the bias voltage as compared to Id and set to a constant. This Poisson’s 
equation is obtained by neglecting the variation of carrier concentration in the direction 
perpendicular to the channel and can be solved with boundary conditions V(L = L c ) — V tat 
and E(L = L c ) = E c . The final form of the solution becomes 


V D - V, at = — - ^ 2 Lc)2 + E c2 (L - L c ) = CI d K l + BKi (12) 

where C = -- = 2f ^Zv,6 third model parameter. Equations (10) and (12) can be 

solved simultaneously to find Ki and V sat : 


K, = 


-X + ^ + [2C^(Pa-P tfe ) 2 -4B][l + ^^]Fc 


CA(V a -V th ) 2 -2B 


(13) 


where X = B+V D + V G -V th . 
be written as 


Then from equation (8), the saturation current equation can 


Id 


MVg ~ Vth - Y r]V 1 at 
1 - K x + y 


(14) 


(B) I-V Equations in the Compressed Transconductance Region 

When the gate voltage is high enough such that V G > V c , the AlGaAs layer starts to conduct 
current. This current conduction mechanism can be considered similar to a parasitic 
MESFET and it is shown in Figure 3. 

When w — 0 at x = L\, from equation (1) the voltage inside the channel is V = V a = 
V p — Vfn + V G , where V p is defined as 




H»o \ 2 




(15) 


When Vp < V 0 , the sheet carrier concentration, n t of the whole 2-DEG channel is 
equal to its equilibrium value, n to and it is independent of gate and drain voltage. The 
2-DEG channel is then like a non-linear resistor with sheet concentration, n to , while the 
undepleted AlGaAs behaves like a MESFET. This equilibrium concentration is assumed 
to be maximum and is given by (from equation (3)) 


n to = (3 {Vk - V p - F tK ) (16) 

iFrom the schematic diagram shown in Figure 3, the compressed transconductance 
region can be divided into three different regions of operation : 

(i) Linear Region I : Vp < V a 



9.36 



i — >• x 


Figure 3: Schematic Diagram Showing Current Saturation in the 2-DEG Channel and the 
Parasitic Conduction through the AlGaAs Layer. 


(ii) Linear Region II : V a < Vp< V fi t 

(iii) Saturation Region : Vp > v m 

Here, the assumption V, a t > V a has been made to allow division into various regions of 
operation. This assumption is true for typical H EM T devices. 

(i) Linear Region I 

For Vp < Vq, the current through the AlGaAs layer is derived as in the case of the MESFET 
and is given by 



T/ 2 (V u - V G + Vpfl 2 - (F w - Vg) 2 ' 2 

v B - - 


(17) 


where E — ( d — di — and F = LE C i are two more model parameters. 

The current through the 2-DEG channel becomes 


r _ A(V bi -V p ~ V th )V D 

1 2 _ - v 


(18) 


The total current in this region of operation is the sum of these two currents : I D = J x + / 3 . 


(ii) Linear Region II 

iFrom Figure 3, for Vp > V„, the current flowing through the 2-DEG channel is 


h = 


ZqfWjoVo 


iFrom this equation 


Li = 


Zqfi 2 n, 0 V 0 V 0 


I 2 E c2 

Current through the AlGaAs layer can be obtained from equation (6) as 


Ii = 


Zq^Njw £ 

~~ n tuv 


(19) 


( 20 ) 
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Integrating this equation for V from 0 to V„ and for x from 0 to L x , and then using equation 
(20) for L\, the final expression for current becomes 


7l ~ i , iH7~ 

h B ~ r fJ v ° 


( 21 ) 


The derivation of the current expression in the 2-DEG is similar to the normal region. 
But here, the limits of integration for V are from V a to Vp and for x from L\ to L. After 
performing the integration and using equation (20) for L x , the current through the 2-DEG 
channel can be obtained as 


, _ A\(V a - V lh - S^)(V D - V.) + (% - V r - V A )V„] 

12 1 j~ Vr. v 22 ) 

1 ~ r B 

So, the total drain current is the sum of equations (21) and (22). 

(iii) Saturation Region 

For the saturation region, Vq > V , at , the current expression for the undepleted AlGaAs 
layer is the same as the linear region II (equation (21)). The principle to find the satu- 
ration voltage in this operating region is similar to that in the normal region except the 
contribution from the parasitic conduction has to be taken into account. From the current 
continuity at the interface of the velocity saturation region and the non-saturation region 
(Figure 3), the saturation voltage can be obtained as 


[(1 - K^BjVg - V th ) + 2 {V G - Vu + V P )V 0 - V?) 
{{V G - V th ) + (1 - K X )B] 


(23) 


On the other hand, the solution of the Poisson’s equation (equation (12)) in this region 
becomes 


V D ~ V sai 


CAK\ [(Vc - V t h)V, at - Sf* - (H, - Hi + H)K + ¥] 

1 -K, + fy 


+ Bh\ 


(24) 


By solving equations (23) and (24) iteratively, values of K x and V, at can be found. Once 
K x and V, al are found, the current through the 2-DEG channel can be obtained as 

t a I(Vg - v lh - *“^)(V,at - Vo) + (V« -V p ~ V th )V 0 ] 

72 - 1-^ + v ( 25 ) 


The total drain current is then the sum of equations (21) and (25). 

In the subthreshold region of operation the charge control is not linear; so, in addition to 
the model and physical parameters, a fitting parameter, D is used to model the threshold 
voltage shift of the 2-DEG caused by the drain voltage. This simple threshold voltage 
correction is given by 

V th = V tho -DxV d (26) 

So, with the nine parameters A, B, C, E, F, V p , V tho) and D, the I-V characteristics 
of the AlGaAs/GaAs HEMT device can be modeled completely. 
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3 Small- Signal Model 


Evaluation and analysis of the small-signal performances of the HEMT are important for 
the operation of microwave circuits. The HEMT is usually biased in the normal transcon- 
ductance region without parasitic conduction for optimal low-noise and/or high-frequency 
performance. Some of the small-signal parameters like transconductance, gate-to-source 
capacitance, current gain cut-off frequency etc. can be derived analytically from this 
model. The derivation of these parameters in the saturated normal region and also in the 
compressed transconductance region is mathematically complicated and computationally 
involves more CPU time. So, to determine these parameters in those regions of operation 
a computationally efficient numerical technique has been used. Methods of determining of 
these small-signal parameters are discussed in the next few subsections. 

3.1 Transconductance, g m 

The intrinsic transconductance, g m at constant drain voltage is defined as 

_ Md , 

9m — Qy I V&— Constant 

The g m in the linear normal region can be obtained analytically by differentiating drain 
current (equation (7)) with respect to gate voltage : 


dVn 


'A(V G -V th -^)V D ' 

l + Yn. 

1 ^ B 


av d 

1 + 

1 ^ B 


(27) 


The transconductance increases with drain voltage before current saturation and is 
inversely proportional to gate length and mobility degradation factor (1 + ^), 

To calculate g m in the saturation region and in the compressed (both linear and sat- 
uration) region, we differentiate the corresponding drain currents numerically. For this 
numerical differentiation we have used the centered-finite-divided difference equation of 
the form [7] 


<7m(F Gi ) 


WG.+J-MVb.,,) 
Va t „ - Vg^ 


(28) 


Here, g m (VGi) is the transconductance evaluated at the i th point. 


3.2 Gate-to-Source Capacitance, C gs 

Gate-to-source capacitance, C g , is defined, with the assumption C g d < C gt , as 

r = 

9 ’ *“ dV G 


where Qt is the total charge. 
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In the normal region, the AlGaAs layer is completely depleted, so the C a . is due only 
to the two-dimensional electron gas. Thus, for the normal region 


' 9 > 


d 

dV G 



Substituting equation (3) for n,(x) and then performing the integration, we get 


C 

9 ‘ dVc 


f q0(V G -V(x)-V th )dx 

Jo 


AL 2 ( 2 + y§-) 


P2 


(29) 


Calculation of gate-to-source capacitance in the saturation region is more complicated 
because of complexity in the total charge calculation in the channel. The method we have 
used to calculate the charge in the channel is given in detail in reference [8]. The final 
expression of total charge, Qt becomes 


Qt 


AL 


(Vb - V tk )L - VUL - y) 


(30) 


where 

, _ V.atL AL(V g - V th - V?)V. at 

c B I D 


(31) 


In this equation, the saturation current, I D is calculated by using equation (14) at the 
saturation voltage, V, at . 

Once we know the total charge in the channel we can calculate the C gt by using nu- 
merical differentiation. The form of this differentiation is analogous to the g m equation 


C g ,(V Gi ) 


Qt(Vg, + 1 ) - Qr(Vb, _,) 

Va w ~ Vb,., 


(32) 


Ideally, to calculate C g3 in the compressed region, the capacitance due to the charge 
accumulated in the undepleted AlGaAs layer has to be added with the capacitance due to 
the 2-DEG channel. But the calculation of the capacitance due to AlGaAs layer analyt- 
ically from this model is not very straightforward. Moreover, this additional capacitance 
contribution may not be very significant, particularly at high drain voltages. So, in this 
work we have neglected this contribution compared to the capacitance due to the 2-DEG 
channel charge. Therefore, equation (32) has also been used to calculate the gate-to-source 
capacitances in the compressed transconductance region. 


3.3 Current Gain Cut-off Frequency, fr 

In microwave applications, the current gain cut-off frequency is the frequency used as an 
indicator of the device speed. The conventional definition of f T is 

f T = JtL- 
J 2tt C g , 


d-r 
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Jn the normal linear region, we calculated fa analytically by using equations (27) and (29): 


1 

AV d 

M 2 

PiV D 

2tt 

1 + Ya. 
L i + bJ 

Uf’(2+T?)J 

- 2*£ 2 (1 + t?)(2 + 50 


(33) 


We again adopted the numerical techniques to calculate fa for norma! saturated region 
and both linear and saturated compressed regions. This numerical expression is given by 


f T (Va) = 

? T[VG ' } 2n C g ,(V G .) 


(34) 


Here, the fo, g m and C g , are calculated at the i th point. 


3,4 Optimum Cut-off Frequency, /t(o pt) 

Another important parameter in microwave applications is the optimum frequency, /r( opt). 
This optimum frequency is defined as the maximum value of the current gain cut-off fre- 
quency just before current saturation occurs. Thus, in the normal transconductance region, 
f T ( opt) is approximated as 


f T {opt) = 


PjV*at 

2ttL2(i + V)( 2 + V) 


(35) 


Here, v. at , the value of the saturation voltage when current just starts to saturate, can be 
evaluated by setting K\ = 0 in equation (10) : 


B(V a ~ Vtk) 

B + (V g - V th ) 


(36) 


4 Results and Discussion 

4.1 The I-V Characteristics 

To validate the dc model we have developed a computer simulation program which cal- 
culates the I-V characteristics over the entire region of operation. Using this simulation 
program we have calculated the I-V characteristics of all the four HEMTs, The device 
physical parameters and the modeling parameters of these HEMTs, taken from reference 
[4], are given in Table 1. 

In the derivation of the drain current equations in section 2, the dc model does not 
include the effects of parasitic source and drain resistances explicitly. These effects can be 
taken into account in the model by solving the nonlinear equations which are given below 

V G s = Vg + Id{Vg,V p )R s (37) 


and 


V ps = V D + Ip{V G ,V D )(Rs + Rd) 


(38) 
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Device 

HEMT #1 (TRW #2078) 

HEMT #2 

HEMT #3 (GE #5410) 

HEMT #4 

L(/«n) 

0.35 

1.0 

0.25 

1.0 

Z(jim) 

65 

145 

100 

1200 

WV) 

-0.017 

-0.901 

-0.912 

-2.389 

V 9 (V) 

- 

- 

1.481 

2.319 

YuiYl. 

- 

- 

0.85 

0.85 

A(mA/ V 2 ) 

49.517 

101.253 

103.539 

454.167 

B(V) 

5.285 

1.604 

0.616 

0.948 

C(K Q) 

8.341 

0.583 

0.992 

0.201 

D 

0.015 

0 

0.092 

0.008 

E(mA/V) 

1 

- 

81.825 

542.663 

F(V) 

' ! 

- 

2.154 

3.04 

Rs{ 0) 

5.9 

7.0 

4.6 

1.0 

R D {n) 

6.0 

7.0 

6.0 

1.0 


Table 1: Physical and Model Parameters of the HEMTs. 


where Vgs and Vds are the externally applied gate and drain voltages respectively; Rs 
and Rd are the parasitic source and drain resistances. These two equations were solved 
iteratively in the program to find the values of Vq and Vp for given values of external 
voltages Vgs and Vds- 

The HEMT #1 and #2 show only normal transconductance effects; only five model 
parameters, V t h 0 , A, B, C and D are needed in the program to calculate the I-V relation. 
With these parameter values and using equations (7), (10), (13), (14), (37) and (38), we 
have developed a simulation program which calculates the drain-to-source current as a 
function of external drain voltage for different external gate voltages. 

Figure 4 shows the I-V curve of the HEMT $1. In the program, we have swept the 
drain voltage from 0 to 3 volts with a 0.2 volts steps and calculated drain-to-source currents 
for gate voltages Vgs — 0, 0.1, 0.2, 0.3, 0.4 and 0.5 volts. As a comparison, we have also 
plotted the experimental data obtained from reference [4]. From the figure, we can see a 
nice agreement between our I-V results and the experimental data. 

Simulated results along with experimental data [4,6] of the HEMT #2 are shown in 
Figure 5. In this case the drain voltage was varied from 0 to 3 volts with 0.25 volts steps. 
Drain currents for Vgs — -0.8, -0.6, -0.4, -0.2 and 0 volts were calculated. The low gate 
bias curves agree very well with the experimental values. As the gate bias increases a small 
deviation occurs near the linear and saturation transition region. 

The I-V characteristics of the HEMT $3 and $4 (double hetero junction HEMT) are 
more complex because of the compressed transconductance effect (in addition to the normal 
transconductance effect). Four additional parameters V p , V^, E and F are needed to model 
this effect. So, with the nine parameter values listed in Table 1 and using the equations 
(17), (18) and (21-25), we have calculated the drain-to-source currents in the compressed 
transconductance region. Equation (24) was rearranged such that K\ can be written in 
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Figure 4: Characteristics of HEMT #1 Figure 5: I-V Characteristics of HEMT #2 
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Figure 6: Characteristics of HEMT #3 Figure 7: TV Characteristics of HEMT $4 
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Device 

Bias Condition j Iif(mA) ] 

9m(mS) | C„(fF) | MGHz) j 

fx{opt)(GHz) 

HEMT #1 

Normal Linear Region 
V GS = 0.4V, V DS = 0.4V 

3.677 

16.526 

28.501 

92,283 

96.055 

HEMT #1 

Normal Saturation Region 
Vas = 0.2V, V DS = 1.0V 

1.376 

10.868 

6.5129 

265.58 

i 

HEMT #2 

Normal Linear Region 
V GS = -0.2V, V DS = 0.5V 

12.389 

27.471 

507.09 

8.622 

10.698 

HEMT #2 

Normal Saturation Region 
V es = -0.2V, V DS = 1.0V 

15.461 

31.098 

134.89 

36.692 

' 

HEMT #3 

Normal Linear Region 
V GS = -0.65V, Vx> s = 0.2V 




|gg||jj| 

69.527 

HEMT #3 

Normal Saturation Region 
V G5 = -0.8V, V DS = 1.0V 

nn 

14.513 



' 

HEMT #3 

Compressed Linear Region 
V GS = -0.2V, V DS = 0.4V 

8.476 

6.699 



- 

HEMT #3 

Compressed Saturation Region 
V GS = -0.4V, V DS = 1.0V 

9.215 

16.924 


ms 

- 

HEMT #4 

Normal Linear Region 
V G5 = -1.5V, V DS = 0.5V 

74.628 

116.277 

1583.0 

11.89 

13.18 

HEMT #4 

Normal Saturation Region 
V GS = -1.5V, V DS = 1.0V 

88.939 

133.306 

500.9 

42.36 


HEMT #4 

Compressed Linear Region 
V GS = -0.5V, V DS = 1.0V 

188.287 

48.913 

1113.0 

6.996 

- 

HEMT #4 

Compressed Saturation Region 
V GS = -1.0V, V 05 = 1.0V 

153.567 

123.322 

548.7 

35.77 



Table 2: Small-Signal Performances of the AlGaAs/GaAs HEMTs Calculated in this Work. 


curves for four HEMTs were successfully calculated and compared with the experimental 
data reported earlier [4,6]. 

In the second phase of the work, analytical and numerical methods were developed 
to predict some of the important small-signal performances of these HEMTs. Based on 
this new computer-aided model, the small-signal parameters, £ m , C gB , fa and fr(opt) 
were calculated and are presented in Table 2. The proposed small-signal model for the 
AlGaAs/GaAs HEMT device may be useful to VLSI and microwave applications in future. 


6 Acknowledgment 

The authors acknowledge NASA for funding this project under the grant NAG5-1043. 
They would like to thank the NASA Space Engineering Research Center, University of 
Idaho for partial funding and for providing facilities to accomplish this work. The second 
author acknowledges Jesus of Nazareth for His sacrifice and example. 





























3rd NASA Symposium on VLSI Design 1991 


9.3.13 


terms of V, at 


_^5 + VD) + V / (^ + VD) 2 + 4[C'^y-B][(l + V)(^D-V.„t)] 

1 2 {CAY - B} 

where Y = (Vg — Vt^V^at — ^ 1 — Equations (23) and (39) were solved iteratively in the 
program by assuming an initial value of K 1 = 0 to obtain V, at and then K x . After K x and 
V, a t are known, the drain current through the 2-DEG channel in the saturation region is 
calculated by using equation (25). 

Figure 6 shows the I-V curve of the HEMT #3. Here, we have scanned the drain 
voltage from 0 to 3 volts at a step of 0.2 volts for the gate voltages, Vbs = -1.0, -0.8, 
-0.6, -0.4, -0.2 and 0 volts. As we can see the modeled result agrees very well with the 
experimental data [4]. 

Finally, we have calculated the I-V characteristics of a double heterojunction HEMT 
(HEMT # 4) and the results, along with the experimental data are shown in Figure 7. 
These results also agree fairly well with the published measured data [4]. 

Two of the four HEMTs (HEMT $1 and $3) are sub-half-micron gate HEMTs. Un- 
modeled short channel effects such as velocity overshoot and unmodeled hot carrier effects 
may occur in these two HEMTs. It is reported that these effects start to become prominent 
below 0.25 pm gate length [9], therefore HEMT ^3 may show considerable short channel 
effect in the compressed transconductance region. Moreover, this dc model was originally 
developed only for the single-heterojunction HEMT, But from our simulation results of 
HEMT #4, which is a double-heterojunction HEMT, we found that this model also appears 
to be good for the double-heterojunction HEMT. 

4.2 Small-Signal Performance Calculation 

Based on the equations derived in section 3 and the physical parameters listed in Table 1, 
we have developed the simulation program which calculates the small-signal performances. 
Using this program we have calculated Id , 5 m , C gt , fo and /r(opt) as a function of gate 
voltage keeping drain voltage fixed. Table 2 shows the small-signal parameter values for 
all the four HEMTs for different drain and gate bias conditions. 

This small-signal model has been developed in an academic environment, based on 
a quasi-static approximation. The values of the small-signal parameters are essentially 
theoretical and have not been rigorously validated in this work because of the unavailability 
of the experimental data. 

5 Conclusion 

A complete analytical dc model for the uniformly doped AlGaAs/GaAs HEMT device 
has extensively analyzed and validated independently. Based on the model a simulation 
program was developed to calculate the I-V characteristics. Using this program, the I-V 
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